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ABSTRACT 

Using a reionization model simultaneously accounting for a number of experimental data 
sets, we investigate the nature and properties of reionization sources. Such model predicts 
that hydrogen reionization starts at z <w 15, is initially driven by metal-free (PopIII) stars, and 
is 90% complete by z w 8. We find that a fraction / 7 > 80% of the ionizing power at z > 7 
comes from haloes of mass M < lO 9 Af0 predominantly harboring PopIII stars; a turnover to 
a PopII-dominated phase occurs shortly after, with this population, residing in M > 1O 9 M 
haloes, yielding / 7 w 60% at z = 6. Using Lyman-break broadband dropout techniques, J- 
band detection of sources contributing to 50% (90%) of the ionizing power at z ~ 7.5 requires 
to reach a magnitude Jho.ab = 31.2(31.7), where ~ 15(30) (PopIII) sources/arcmin 2 are 
predicted. We conclude that z > 7 sources tentatively identified in broadband surveys are 
relatively massive (M « 1O 9 M ) and rare objects which are only marginally (« 1%) adding 
to the reionization photon budget. 
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1 INTRODUCTION 

The study of reionization received a big boost due to availability 
of a variety of observational data accumulated over the past few 
years (for reviews see Fan, Carilli, & Keating 2006 and Choudhury 
& Ferrara 2006a); additional progresses are soon expected from a 
number of different ground-based (LOFAR, MWA, PAST, SKA, 
ALMA) and space-born (JWST, PLANCK, GLAST) experiments. 

The available results have allowed us to build self-consistent 
reionization scenarios that are able to account simultaneously for a 
number of observables (redshift evolution of Lyman-limit absorp- 
tion systems, GP and electron scattering optical depths, mean tem- 
perature of the intergalactic medium (IGM), and cosmic star forma- 
tion history; Choudhury & Ferrara 2005, 2006; hereafter CF05 and 
CF06 respectively). To summarize the emerging picture, the most 
favorable model is one in which hydrogen reionization was an ex- 
tended process starting around z « 15 and being 90% complete by 

2 « 8. This (early) reionization model was also shown not to be 
in conflict with the Gunn-Peterson optical depth evolution deduced 
from QSO absorption line experiments at z > 6 both by using 
the statistics of dark gaps in the Lya transmitted flux (Gallerani, 
Choudhury, & Ferrara 2006) and through radiative transfer simu- 
lations of ionized regions around QSOs (Maselli et al. 2007). Ac- 
cording to the model, reionization is initially driven by metal-free 
stars in low mass (M < 10 s Mq) haloes; the conditions for the for- 
mation of these objects are soon erased by the combined action of 
chemical and radiative feedback at z < 10. As a consequence, the 
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photoionizing power (and therefore integrated luminosity) of these 
sources is most significant around z ~ 8 — 12. 

In spite of this successful overall picture, relatively less at- 
tention has been devoted so far to the observable properties of the 
sources responsible for cosmic reionization. Our main aim in this 
work is to fill such gap by providing quantitative guidelines for 
observers searching for the first cosmic light sources. Specifically, 
we use the CF05 and CF06 model to estimate the IR fluxes and 
magnitude-limited counts of the primary reionization sources. 

2 BASIC FEATURES OF THE MODEL 

The main features of the semi-analytical model used in this work 1 
could be summarized along the following points (for a more de- 
tailed description see CF05 and CF06). The model accounts for 
IGM inhomogeneities by adopting a lognormal distribution accord- 
ing to the method outlined in Miralda-Escude, Haehnelt, & Rees 
(2000); reionization is said to be complete once all the low-density 
regions (say, with overdensities A < A cr i t ~ 60) are ionized. The 
mean free path of photons is thus determined essentially by the 
distribution of high density regions. We follow the ionization and 
thermal histories of neutral, HII and Helll regions simultaneously 
and self-consistently, treating the IGM as a multi-phase medium. 
Three types of reionization sources have been assumed: (i) 

1 Throughout the paper, we use the best-fit cosmological parameters from 
the 3-year WMAP data (Spergel et al. 2006), i.e., a flat universe with Q m = 
0.24, n A = 0.76, and U b h 2 = 0.022, and h = 0.73. The parameters 
defining the linear dark matter power spectrum are erg = 0.74, n s = 0.95, 
dn s /dlnfe = 0. 
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Figure 1. Left: Number of ionizing photons per H-atom contributed by haloes of different mass in a fraction of the Hubble time tu (z) equal to the recombi- 
nation time i rec (2) as a function of redshift z for different halo masses and for different stellar populations. The solid and dot-dashed pair of curves, from top 
to bottom, represent haloes of masses 10 7 M Q < M < 1O 8 M , 10 8 M Q < M < 1O 9 M , and M > 10 9 M Q respectively. The dashed line represent the 
evolution of the volume filling factor Qhii of ionized regions. Right: Cumulative fraction of the ionizing power / 7 contributed by haloes of mass > M. The 
curves from right to left correspond to 2 = 6, 7, 8, 9, 10 respectively. 



metal-free (i.e. PopIII) stars having a Salpeter IMF in the mass 
range 1 — IOOMq: they dominate the photoionization rate at high 
redshifts; (ii) PopII stars with sub-solar metallicities also having a 
Salpeter IMF in the mass range 1 — 100M Q ; (iii) QSOs, which are 
significant sources of hard photons at 2 < 6; they have negligible 
effects on the IGM at higher redshifts. Note that there is no com- 
pelling reason to rule out the possibility that the IMF of PopIII stars 
was top-heavy; however, as this is not required by current data, we 
limit ourselves to the most conservative model. 

Reionization by UV sources is accompanied by photo-heating 
of the gas, which can result in a suppression of star formation 
in low-mass haloes. We compute such (radiative) feedback self- 
consistently from the evolution of the thermal properties of the 
IGM. Furthermore, the chemical feedback inducing the PopIII — ► 
PopII transition is implemented according to the detailed study by 
Schneider et al. (2006) in which a merger-tree "genetic" approach 
is used to determine the termination of PopIII star formation in a 
metal-enriched halo. 

The predictions of the model are compared with a wide range 
of observational data sets, namely, (i) redshift evolution of Lyman- 
limit absorption systems (Storrie-Lombardi et al. 1994), (ii) IGM 
Lya and Ly/3 optical depths (Songaila 2004), (iii) electron scat- 
tering optical depth (Spergel et al. 2006), (iv) temperature of the 
mean intergalactic gas (Schaye et al. 1999), (v) cosmic star forma- 
tion history (Nagamine et al. 2004), and (vi) source number counts 
at 2 « 10 from NICMOS HUDF (Bouwens et al. 2005). 

The best-fit reionization model is characterized by a PopII 
(PopIII) star-forming efficiency e*,n = 0.1(e*,m = 0.03) and 
escape fraction / eB c,ii = 0.01(/ eB c,iii = 0.68) (keeping in mind 
that /esc, 11 and / eB c,in are not independent). 2 The resultant value 
of the electron scattering optical depth is T e i = 0.1 

The data constrain the reionization scenario quite tightly. We 



2 The values are slightly different from those quoted in CF06 because of an 
improved likelihood analysis. The qualitative results and main conclusions 



find that hydrogen reionization starts at 2 ~ 15 driven by metal- 
free (PopIII) stars, and it is 90 per cent complete by 2 ~ 8. This 
can be seen from the dashed curve in Fig. 1 which represents the 
evolution of the volume filling factor Qhii (2) of ionized regions. 
After a rapid initial phase, the growth of the volume filled by ion- 
ized regions slows down at 2 < 10 due to the combined action of 
chemical and radiative feedback, making reionization a consider- 
ably extended process completing only at 2 « 6. 

3 PROPERTIES OF THE REIONIZATION SOURCES 

Having identified the most probable reionization scenario, we 
can now confidently determine the properties of the reionization 
sources governing it. Let us start by defining the quantity 
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as the number of ionizing photons per H-atom contributed by 
haloes in the mass range [M m i n , M max ] in a fraction of the Hubble 
time tniz) equal to the recombination time t rcc (z). By construc- 
tion, the IGM is reionized when x-, > 1. The quantity uh is the 
comoving number density of hydrogen atoms while n 7 is the time- 
integrated comoving photon density, calculated using the relation 
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are unaffected by this modification. 



where n 7 (M m i n : M max , t) is the ionizing photon comoving emis- 
sivity from haloes within [M m i n , M max ] . 

The plot of Xj(z) for different mass ranges is shown in Fig. 
1. The lower mass 1O 7 -1O 8 M0 haloes dominate the photon pro- 
duction rate at early redshifts providing about 0.25 photon/H-atom 
on the fractional recombination timescale. These objects produce 
the first ionized regions, are preferentially metal-free, and there- 
fore mostly harbor PopIII stars of high specific ionizing power. 
At 2 < 8, the contribution from PopIII haloes decreases because 
their formation is hampered by the heating associated with radia- 
tive feedback. As a result, the progress of ionization fronts relies on 
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7<z<8, J no Band 9.5<z<10.5, H 160 Band 




Figure 2. Number density of reionization sources (PopII and PopIII) as a function of the limiting magnitude, m\ B , at 7 < z < 8 observed in Jno band (left 
panel) and at 9.5 < z < 10.5 observed in Hieo band (right). Halo masses corresponding to tuab for the particular redshift under consideration [given by 
equations (4) and (5)] are shown by the straight lines and the relevant values can be read off from the right vertical axis. 



photons emitted by more massive haloes with M > 10 9 M Q . It can 
be seen from Fig.l that such high mass haloes do not host PopIII 
stars as they form from the merging of already polluted progeni- 
tors, a result of the "genetic" transmission of chemical feedback 
(Schneider et al. 2006). This combination of radiative and chem- 
ical feedback makes the reionization process quite extended and 
its completion has to wait until z « 6 when the PopII stars (and 
QSOs, not shown in the Figure; see CF06 for details) dominate the 
photoionization rate. 

Additional insights on the source properties may be gained 
by considering the fractional instantaneous contribution of haloes 
above a certain mass, 



n 7 (> M, z) 
n 7 (z) 



(3) 



also shown for different redshifts in the right hand panel of Figure 
1 . The results shown emphasize again that > 80% of the ionizing 
power at z > 7 is provided by haloes with masses < 1O 9 M0 
which are predominantly harboring PopIII stars. A turnover to a 
PopII-dominated reionization phase occurs shortly after, with this 
population, residing in M > 10 9 M© haloes, producing « 60% of 
the ionizing photons at z — 6. In conclusion, PopIII stars and small 
galaxies initiate reionization at high redshift and remain important 
until they are overcome by PopII stars and QSOs below z = 7. 

4 SOURCE COUNTS AT HIGH REDSHIFT 

Having determined which haloes contribute most significantly to 
the ionizing power at a particular redshift, we now address the de- 
tectability of these sources in broadband imaging surveys and the 
optimal strategies to do so. A halo of (dark matter) mass M at a 
redshift z will emit a flux given by (Salvaterra & Ferrara 2006) 



4nd 2 L (z')Avo 



J dv l v ,(M) • 



-T e ff (fo,z=0,z') 



(4) 



where e, is the star-forming efficiency of the population under con- 
sideration, l„i (At) is a template specific luminosity for the stellar 



population of age At = t x i —t z it (the time elapsed between the two 
redshifts), cLl{z') is the luminosity distance and Ai/rj is the instru- 
mental bandwidth. The quantity T ct t(vo, z = 0, z') is the effective 
optical depth due to intervening gas at vq between z' and z = 0. 
The main contribution to T e g (vo, z = 0, z') comes from the com- 
bined blanketing of the Lyman series lines in the emitter rest frame 
wavelength range 912A < A < 1216A and from the continuum 
absorption from neutral hydrogen in the range A < 912A. Both the 
contributions to the opacity can be calculated self-consistently from 
the (modified) lognormal density distribution (CF05) assumed in 
our model. For calculating the template luminosity l v , we use stel- 
lar population models (having metallicity Z = 0.004 = 0.2Z & ) of 
Bruzual & Chariot (2003) for PopII stars and of Schaerer (2002) for 
PopIII stars. Note that the star-forming efficiency e* in the above 
equation is fixed by constraining the reionization history and is not 
a free parameter as far as the calculation of the source counts in this 
work is concerned. The flux F„ can be transformed into a magni- 
tude in the AB system: 



m A B = -2.5 log! 
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-48.6 



(5) 



The number of sources within a redshift interval [z m i n , z m ax] 
observed in a solid angle dfi having a flux larger than F VQ is 



(6) 



where dV/dz'dQ denotes the comoving volume element per unit 
redshift per unit solid angle, and 



dFl- {F ^ Z)= J z , dZ dK n iF ^ At) dMd^ (M ' Z ] 



(7) 



is the comoving number of objects at redshift z' with observed flux 
within [F' vo , F' VQ + dFl Q ] . The quantity d 2 n/dMdz" gives the for- 
mation rate of haloes of mass M, which we obtain from our reion- 
ization model. 
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For definiteness, let us consider two bands corresponding to 
the NICMOS observations of the HUDF (Bouwens et al. 2004; 
Bouwens et al. 2005), namely Jno and ffieo; these broadband 
filters are appropriate to detect Lyman-break dropout sources at 
z ~ 7.5 and z ~ 10, respectively. The main results of the cal- 
culation are shown in Fig. 2 where we plot the number of sources 
(i.e. galaxies powered by either PopII or PopIII stars) observed in 
a particular redshift interval and band as a function of the limiting 
magnitude tuab- The halo masses corresponding to m\B for the 
particular redshift under consideration [given by equations (4) and 
(5)] are shown by the straight lines and the relevant values can be 
read off from the right vertical axis. 

From Fig. 2, we deduce that, at the currently achieved sen- 
sitivity limit of ~ 28 AB magnitude, our best-fit model predicts 
~ 10(1) Lyman-break sources per arcmin at z ~ 7.5(10) observ- 
able in the Jno (-H160) band. Note that all these sources are bright 
PopII star forming haloes having masses ~ 10 10 M Q . In the previ- 
ous Section we found that such sources provide only a negligible 
(~ 1%) contribution to reionization at 2 > 6. 

Some of these sources have been tentatively identified in 
broadband observations. For example, observations of HUDF using 
the NICMOS filter (Bouwens et al. 2004; Bouwens et al. 2005) at 
AB magnitude limit ~ 28 reveal ~ 1 source per arcmin 2 at z ~ 7.5 
and < 1 source per arcmin 2 at 2 ~ 10 respectively. On the other 
hand, deep near-IR photometry of two lensing clusters (A 1835 and 
AC 114) obtained with ISAAC/VLT (Richard et al. 2006) seem to 
indicate ~ 1 source per arcmin 2 at 2 ~ 7— 10 at a magnitude 
limit equivalent to _Hi6o,ab ~ 26, which is considerably higher 
than the NICMOS results. Stringent constraints are difficult to ob- 
tain from theoretical models as there remains considerable confu- 
sion regarding the actual number of reliable detections. However, 
although the actual source count is disputed, it is most likely that 
these tentatively detected sources are bright massive haloes which 
are not significant for reionization. 

The sources responsible for the bulk of reionization are galax- 
ies formed inside low mass < 10 M© haloes and powered by 
PopIII stars. By inspecting the right hand panel of Fig. 1 we con- 
clude that in order to observe sources which contribute to 50% 
(90%) of the ionizing power at 2 ~ 7.5, it is necessary to observe 
PopIII sources of mass < 1O 81 M (< 10 7 ' 9 M Q ). This in turn re- 
quires to reach a magnitude sensitivity of Jho,ab = 31.2(31.7) 
[left panel of Fig. 2], currently beyond the reach of any instru- 
ment. Interestingly, one can expect to observe ~ 15(30) PopIII 
sources per arcmin 2 once such sensitivities are reached. Similar 
conclusions can be drawn from the right hand panel for sources at 
2 ~ 10 observable in the ffieo band. For detecting sources which 
contribute to 50% (90%) of the ionizing power, the magnitude sen- 
sitivity required would be -Hi60,ab = 32.2(32.5), which would 
again correspond to halo masses of 1O 7 ' 8 M Q (1O 7 ' 7 M ). Hence, 
the detection of reionization sources would require sensitivities bet- 
ter than 32 AB magnitudes in this band. 

4.1 Variants of the best-fit model 

In this subsection, we try to obtain some indication about how 
much can the source count estimates vary without violating any 
other observational constraints. We first note that the star-forming 
efficiency of PopII sources are quite stringently constrained by low- 
redshift observations, e.g., the cosmic SFR and the QSO absorption 
lines; hence their numbers cannot vary significantly. The situation 
is markedly different for PopIII stars as there are very few observa- 
tions constraining their properties. In order to calculate the varia- 
tion in PopIII source counts, we consider two extreme limits of the 
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Figure 3. Number density of PopIII sources as a function of the limiting 
magnitude, mAB, at 9.5 < z < 10.5 observed in H16O band. The plots 
are for three models of PopIII star-formation as discussed in the text. Halo 
masses corresponding to m^B for the particular redshift under considera- 
tion [given by equations (4) and (5)] are shown by the straight lines and the 
relevant values can be read off from the right vertical axis. 

parameter e»,m which are essentially determined by the bounds on 

(i) The high-e»,in model: This model is characterized by the 
parameters e»,n = 0.1, e»,m = 0.15, / eS c,ii = 0.0, / csc ,iii = 0.4. 
The resultant r e i is 0.12 (the \-a upper limit given by WMAP3; 
Spergel et al. 2006). The initial stages of reionization in this case 
proceeds much faster than in the best-fit model and 90 per cent of 
the IGM is ionized by 2 w 10.5. However, because of the high 
efficiency of the PopIII stars, the radiative feedback is more severe 
and hence the reionization gets extended till 2 « 6 as in the best-fit 
model. 

(ii) The low-e*,m model: This model is characterized by the 
parameters e*,n = 0.1,e.,m = 0.01, / csc ,n = 0.04, / csc ,iii = 
0.06. The resultant r e i is 0.06 (the 1-er lower limit given by 
WMAP3; Spergel et al. 2006). This model corresponds to the min- 
imum ionizing power contribution required from PopIII stars. This 
model has an important qualitative difference from the high-e*,m 
model; because of such low PopIII efficiency, almost all of the 
reionization process is driven by the PopII stars. In fact, the lack 
of ionizing power at high-2 makes reionization start much later. 
Also, the effects of feedback are minor in this model; as a result, 
the reionization occurs rapidly and shortly before 2 = 6. 

We now proceed to discuss how these models differ in their 
predictions for PopIII source counts (Fig. 3). Since the effects of 
PopIII stars are most noticeable at 2 > 10, we focus on the Hieo 
band. The first important point is that the number of bright sources 
increases with e»,m, a direct consequence of the higher specific 
SFR within dark matter haloes. This would imply that one can 
observe most of the reionization sources at a magnitude limit of 
#160, ab ~ 30.5, a limit almost certainly attained by near future 
experiments. However, the number of sources detected at fainter 
magnitudes is smaller for the high-e*,in model than for the other 
two (best-fit and high-e*,in) models. This is because the efficient 
star formation at high-2 makes the radiative feedback more effec- 
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tive, hence effectively quenching star formation in low-mass haloes 
and suppressing the number of faint sources. Of course, one finds 
a rise in the number of faint sources for the low-e*,m model when 
the feedback has the least effect; one should be able to observe, in 
principle, > 100 sources at -ffieo.AB ~ 33. However, one should 
keep in mind that for low values of e*,m w 0.01, PopIII stars are 
not the main driving force for reionization anyway, and the main fo- 
cus should be the source counts for PopII sources. In fact, we find 
that sources contributing to 90% of the ionizing power at z £s 7.5 
would be the PopII stars with haloes of mass < 10 s Mq\ they can 
be observed with a sensitivity limit of Jho.ab ~ 33, which is quite 
difficult with near-future experiments. 

It is thus clear from the above discussion that the luminos- 
ity function of PopIII sources at 2 ~ 10 could be important for 
constraining the star-forming efficiency of PopIII stars and also for 
some indirect understanding of feedback effects. 

5 DISCUSSION 

We have used the self-consistent model of CF06, which is con- 
sistent with a variety of observations, to estimate the number of 
sources (galaxies) at 2 ~ 7 — 10. According to our analysis: 

• The best-fit model predicts > 1 sources at 2 ~ 7 — 10 per 
arcmin 2 at a sensitivity limit of ~ 28 AB magnitude, roughly con- 
sistent with present observations. However, these sources are bright 
massive (> 1O 1O M0) haloes forming PopII stars which have neg- 
ligible contribution to reionization of the IGM. 

• Reionization at 2 > 6 is actually driven by PopIII stars in low 
mass (< 10 8 M Q ) haloes. The required sensitivity to detect these 
sources would be 31-32 AB magnitude. 

• In case the star-forming efficiency of PopIII stars is higher 
than what is assumed in our best-fit model, the reionization sources 
can be detected with sensitivities of 30 AB magnitude; however, 
the number of sources detected could be smaller than the best-fit 
model because of negligible star formation in low-mass haloes due 
to radiative feedback. 

• For low star-forming efficiencies of PopIII stars (~ 0.01), the 
reionization is actually driven by PopII stars and occurs rapidly and 
shortly before 2 = 6. In this case, the sensitivity required to ob- 
serve the reionization sources (PopII stars within haloes of masses 
< 10 s M ) at 2 > 6 would be 33 AB magnitude. 

The sensitivities required to observe the reionization sources 
are expected to be achieved in future deep imaging surveys, partic- 
ularly the ones which take advantage of gravitational lensing mag- 
nification. For example, the NIR Wide Field Camera 3 (WFC3), 
scheduled to be installed on HST in near future, promises to achieve 
a sensitivity limit of mAB ~ 31 in the Jno and I/160 filters for a 
field of view of ~ 5 arcmin 2 in about a few hundred hours of obser- 
vation time (Stark, Loeb, & Ellis 2007). According to our estimates, 
such lensed surveys should be able to detect quite a few PopIII 
reionization sources (via Lyman-break dropout techniques) in the 
field of view. A much better prospect of detecting these sources 
would be through the Ultra-Deep Imaging Survey using the JWST 
which too plans to achieve a sensitivity limit of ttiab > 31 over 
100-200 hours of observation time per filter (Gardner et al. 2006). 
Direct detection of these sources would put stringent constraints on 
reionization history, and in addition can be used for understanding 
physical processes like feedback. 

Provided that the contamination problems due to bright atmo- 
spheric emission lines, which could restrict visibility up to 50% of 
the redshift range in the J-band, could be maintained under control, 



a complementary approach for detecting the reionization sources 
would be through narrow band surveys for Lya emitters. A rough 
estimate of the Lya luminosity from a halo of mass M forming 
PopIII stars is given by 

L a = e*, m(l - j 'esc, in) ^p- M c-Ly a q{H) (8) 

where CL yQ1 = 1.04 x 10 11 ergs and q(H) is the rate of hydrogen- 
ionizing photons per unit mass of stars formed (Schaerer 2002). 
The above relation does not take into account the attenuation aris- 
ing from neutral hydrogen around the source, hence the luminosi- 
ties could possibly be overestimated. Under the above assump- 
tions, a 10 8 M Q halo would produce a luminosity of ~ 10 41 
erg s _1 in our best- fit model. Such luminosities seem to be well 
within the reach of lensed Lya surveys. For example, Stark et al. 
(2007) have detected ~ 2 sources having luminosities ~ 10 41 ' 5 
erg s _1 at 2 ~ 10 within a area of 0.3 arcmin 2 using a Keck 
survey of gravitationally-lensed sources. There is a possibility that 
these sources could be PopIII stars forming in haloes of masses 
> 10 8 Mq which do have < 30% contribution to the ionizing 
power at 2 « 10. Future narrowband Lya surveys like DAzLE 
thus would be extremely important for direct detection of primary 
reionization sources at high redshift. 

Finally, it is worth mentioning that our models do not include 
(i) star formation in minihaloes where molecular cooling is effi- 
cient, or (ii) the possibility of a top-heavy IMF for PopIII stars. The 
predictions regarding the source counts could vary considerably de- 
pending on the details of the above two processes, and hence future 
surveys present an excellent opportunity to probe such effects. 
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